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Целью работы является описание теоретических и экспери- 
ментальных исследований коэффициента ослабления вибра- 
ций от колесных пар на основание кузова электровоза. Изуче- 
ние виброакустических характеристик на рабочих местах ло- 
комотивных бригад показало, что превышения октавных 
уровней звукового давления над предельно допустимыми в 
значительной степени определяются структурной составляю- 
щей шума, которая формируется вибрациями элементов ку- 
зовных конструкций. Расчет скоростей колебаний элементов 
кузова основан на энергетических методах, традиционно при- 
меняемых для энергетически замкнутых систем. Система 
уравнений учитывает не только геометрические и физико- 
механические параметры элементов, но и вводимую в основа- 
ние кузова вибрационную мощность от колесных пар при 
движении локомотивов. Для инженерных расчетов спектров 
вибрации элементов кузова и, соответственно, уровней звуко- 
вого давления при оценке вводимой вибрационной мощности 
целесообразно использовать частотнозависимую величину 
коэффициента ослабления вибрации. Это особенно важно 
учитывать на этапе проектирования машины. Приведены ре- 
зультаты экспериментальных исследований указанной вели- 
чины и вывод на этой основе регрессионных зависимостей для 
различных скоростей движения. 


Ключевые слова: вибрация, колесная пара, элементы кузова, 
электровоз, коэффициент ослабления вибрации. 
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Введение. Изучению закономерностей формирования спектров вибрации и шума при движении подвижного 
состава посвящены работы [1-11]. Теоретические и экспериментальные исследования виброакустических характери- 
стик на рабочих местах локомотивных бригад, приведенных в работах [10, 11], показали, что превышение предельно 
допустимых значений октавных уровней звукового давления в значительной степени определяется структурной со- 
ставляющей шума. 

Обработка результатов экспериментальных исследований. Для уточнения расчета виброскоростей и 
уровней структурного шума проводились экспериментальные исследования коэффициентов ослабления вибраций 
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(Кос), передаваемых от колесных пар на пол кузова. Методика проведения экспериментов подробно описана в работе 
[10]. При этом рассматривался более широкий низкочастотный диапазон, поскольку нормирование вибрации на рабо- 
чем месте машинистов проводится в этой области частот. При проведении экспериментов выяснилось, что коэффици- 
ент ослабления вибрации зависит не только от их частоты, но и от скорости движения электровоза. 

Исследования вибрации представляли собой серию экспериментов, в ходе которых значения коэффициента 
ослабления четырежды измерялись согласно стандартной матрице планирования двухуровневого двухфакторного 
ПФЭ (полный факторный эксперимент) 2? (табл. 1). В последнем столбце табл. 1 указан порядок (на основе генерато- 
ра случайных чисел) выполнения отдельных опытов при их рандомизации по времени. 






































Таблица 1 
План ПФЭ 2? 
№ Факторы 
Порядок реализации опытов 
п/п х Хх 
1 = - 4 
2 + 2 
3 в + 1 
4 + + 3 
Результаты измерений коэффициента ослабления вибрации представлены в табл. 2. 
Таблица 2 


Результаты экспериментов 

















о — — = = п — 2 
— р Гц | Аа | Же | а | Я, | (и-Е,) | (&2-Е,) | (к-Е,) Х (вы вы) 
1| 70| 8 [0,126 [0,132 | 0,120 | 0,126 | 0,000Е+00 | 3,974Е-05 | 3,974Е-05 7,949Е-05 
2 [10| 8 | 0,204 0,215 | 0,194 | 0,204 | 0,000Е+00 | 1,044Е-04 | 1,044Е-04 2.088Е-04 
3 | 70 [8000 | 0,180 [0,189 | 0,171 | 0,180 | 7,7048-34 | 8,139Е-05 | 8,1398-05 1,628Е-04 
4 | 100 [8000 | 0,161 | 0,169 | 0,153 | 0,161 | 7,704Е-34 | 6,470Е-05 | 6,470Е-05 1,294Е-04 









































Для обработки результатов экспериментов необходимо определить математическое ожидание и среднеквад- 
ратичное отклонение по формулам: 


р == ы Е . (1) 
и=1 
п — 2 
в = У [ны -й. . (2) 
1 


и= 


Проведем оценку статистической однородности результатов эксперимента по критерию Кохрена: 


Ш —_ 2 
х (®.. ви — Кыа 





и= тах 
= < 
р №М тп пи 2 - 6. . (3) 
у у [в =!) 
= и=1 
Здесь { — номера серий опытов (от 1 до №); и — количество параллельных опытов в серии (от | до п); С» и С„ — рас- 
четный и табличный критерии Кохрена для 5-процентной значимости при М сравниваемых и п параллельных опытах: 
0,00021 
р==——=0,36<(,=0,76, 
0,00058 


Таким образом, считаем результаты статистически однородными, те. соответствующими нормальному рас- 
пределению. 
Определим дисперсию воспроизводимости В натуральных переменных по формуле: 
ЕЕ) 
ти 25 себ _ 0,00058 
№("-1) 4(3—1) 





=7,26.105. (4) 


Результаты измерения коэффициентов ослабления вибрации (их средние значения) представлены на рис. 1. 


Подуст С. Ф. Оценка коэффициента ослабления вибраций от колесных пар на основание кузова электровоза 
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Рис. 1. Коэффициенты ослабления вибраций при скоростях движения: 1 — 80 км/ч; 2 — 70 км/ч; 
3 — 90 км/ч; 4 — 100 км/ч 


Определение регрессионных зависимостей. Анализ результатов экспериментов показывает, что максималь- 
ные значения коэффициента ослабления вибраций находятся в среднечастотной области и имеют четкую зависимость 
от скорости движения локомотива. При этом просматривается интересная особенность — с увеличением скорости 
движения локомотива экстремум кривых смещается в область более низких частот. Так, для скорости 70 км/ч экстре- 
мум находится на частоте 2000 Гц, а при скорости 100 км/ч — 125 Гц. Эта особенность осложняет поиск аналитиче- 
ской зависимости функции частотнозависимого коэффициента ослабления вибрации, поэтому принято приближение 
экстремума всех кривых к частоте 500 Гц, поскольку разница между экспериментальными и расчетными данными 
незначительна. 

Изменение коэффициента ослабления вибрации в зависимости от скорости движения также позволяет вы- 
явить экстремальное значение, особенно в среднечастотной области, на скорости 80 км/ч. Таким образом, функция 
экстремальна по двум координатам, и ее можно аппроксимировать как функцию двух переменных, так как выполня- 
ется условие ортогональности, обеспечивающее получение независимых оценок коэффициентов регрессии и их дис- 
персий [12]. 

Экстремальность функции в обеих плоскостях требует применения квадратичной аппроксимации или нахож- 
дения степенной функции адекватного поведения. При этом необходимо учитывать, что один из факторов, а именно 
частота }, является величиной геометрического ряда, последующие значения которого увеличены в 2 раза по сравне- 
нию с предыдущим. Для смягчения влияния этого параметра его необходимо прологарифмировать. 

Таким образом, учитывая экстремальность функции и необходимость логарифмирования входных факторов, 
принято решение аппроксимировать зависимость степенными функциями. Для этого можно использовать метод пере- 
вала, который состоит в построении двух полнофакторных экспериментов таким образом, что значение факторов, 
обеспечивающих оптимум параметра выхода, является общей границей интервалов их варьирования, поставленных с 
двух сторон оптимума. 

Границей интервалов являются значения факторов, находящихся в экстремуме функции. В данном случае 
определены следующие интервалы: 

— для первой области — {= 8—500 Гц, У = 70—80 км/ч; 
— для второй области — {= 500-8000 Гц, У = 80-100 км/ч. 

Для крайних точек интервалов необходимо определить выход, т. е. значение коэффициентов ослабления виб- 
рации. Они комбинируются следующим образом (табл. 3): 














Таблица 3 
Комбинация выходов при ПФЭ 2? 
м №2 
шт шт 
тах шт 
шт тах 
тах тах 
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В нашем случае такая комбинация имеет следующий вид (табл. 4): 
































Таблица 4 
Комбинация выходов для обеих областей 
У, км/ч Гц Кос У, км/ч Гц К 
70 8 0,126 80 500 0,311 
80 8 0,176 100 500 0,257 
70 500 0,25 80 8000 0,215 
80 500 0,311 100 8000 0,161 
































Для получения степенной зависимости необходимо прологарифмировать все параметры модели. Учитывая 
порядок коэффициента ослабления вибрации, искусственно увеличим его на два порядка. Тогда исходные данные бу- 
дут следующими (табл. 5). 






























































Таблица 5 
Исходные данные для расчета коэффициентов регрессии 
ШУ пу 10(100-А„.) ШУ ну 1(100-К,.) 
4,248495 2,079442 2,533697 4,382027 6,214608 3,437208 
4,382027 2,079442 2,867899 4,60517 6,214608 3,246491 
4,248495 6,214608 3,218876 4,382027 8,987197 3.068053 
4,382027 6,214608 3,437208 4,60517 8,987197 2,778819 
При этом уравнение регрессии в кодовых переменных будет иметь вид: 
у=В-Ьлх, +Ь.х. +6 хх, (5) 


где у=1(100-^„.); х=МУ; х. =ШМХ. 
Для нахождения коэффициентов уравнения регрессии необходимо скомпоновать выходы так, как представле- 
но в табл. 6. 

















Таблица 6 
Компоновка выходов для расчета коэффициентов регрессии 
У —У —У У 
У У у —У 
У У У —У 
У У У У 
Хо р 52 52 























У У у у 
Тогда № ==; БЕ; ==: р. =. 
да 9% 4 ОИ а" 274 


Таким образом, на основе выше полученных данных, уравнения регрессии в кодовых переменных имеют вид: 
у = 3,014 + 0,138, +0,314х, — 0,029хх.; у=3,133-—0,12х, —0,209х. —0,025х,х>. 

Определим значимость коэффициентов уравнений регрессии сравнением их абсолютного значения с довери- 
тельным интервалом: 


А-а, - 92 < В]. (6) 
Здесь , — критерий Стьюдента для числа степеней свободы у = № (п-1) (в данном случае приу=8и 1, =2,3 ); 


2 в 
$, — дисперсия коэффициентов уравнений регрессии: 


5%. 7.9610” 


= 6,05.10- 8. 
№-п 4.3 





52 = 


Тогда ДЬ = 2,3-/6,05-107 = 0,057. 


Анализ полученных коэффициентов уравнений регрессии приведен в табл. 7. 


Подуст С. Ф. Оценка коэффициента ослабления вибраций от колесных пар на основание кузова электровоза 












































Таблица 7 
Значимость коэффициентов уравнения регрессии 
ДЬ |, | Вывод ДЬ |, | Вывод 
0 < 3,014 Значим < 3,133 Значим 
1 < 0,138 Значим < 0,120 Значим 
7 7 
р) 0.05 < 0,314 Значим 0.08 < 0,209 Значим 
12 > 0,029 Незначим > 0,025 Незначим 









































Из табл. 7 видно, что коэффициенты при произведении х! и х› статистически незначимы, поэтому данными 
членами уравнений можно пренебречь. Тогда 
у=3,014+0,138ж + 0,314х. ; у = 3,133- 0,12х, —0,209х». 


Для перехода от кодовых переменных к натуральным необходимо выполнить следующие замены перемен- 
НЫХ: 


2(МУ МУ} +1. 





з 


м 





МУ ах - МУ 
2(ШЕ- . 
= (пу п нах) | 


т нии —т т 


Тогда для первого уравнения: 
2(тУ —4,382) 























х = +1 = 14,978 тУ — 64,63, 
4,382 —4,2485 
2(ту-6,215) 
х = +1 = 0,484 т } —2,006. 
6,215—2,079 
Для второго уравнения 
2(тУ -—4,605) 
м = +1 = 8,963 У — 40,275 , 
4,605 —4,382, 
2(ту-8,987) 
хо = +1= 0,7211 } —5,483. 


8,987 —6,215 
Подставляя замену и приводя подобные, получаем уравнения с натуральными переменными: 
т(100-К„)=-6,543+2,069тУ + 0,1521 Г ;т(100-К„,)=9,112-—1,075тУ —0,151т { . 
Проводим потенцирование и получаем степенные зависимости 
к. =Ь4-10-5у2989 =0152 , к. =90, 657-1075 ;-б1я 
Для оценки погрешности полученных зависимостей сравним расчетные и экспериментальные значения коэф- 
фициентов ослабления вибрации (табл. 8). 

















Таблица 8 
Погрешность расчетных и экспериментальных значений 
У, км/ч ГГц Кое Коер А, % У, км/ч ГГц Го Кос А, % 
70 8 0,126 0,129 2,94 80 500 0,311 0,319 2,49 
80 8 0,176 0,171 —2,86 100 500 0,257 0,251 —2,43 
70 500 0,25 0,243 —2,86 80 8000 0,215 0,210 —2,43 
80 500 0,311 0,32, 2,94 100 8000 0,161 0,165 2,49 









































Наибольшая погрешность полученных зависимостей не превышает 3 %. 

Графическая интерпретация этой модели представляет собой две поверхности, соединенные в общей точке 
экстремума исследуемой функции. Для построения этих поверхностей необходимо определить значения промежуточ- 
ных точек по полученным зависимостям. Точка экстремума не является центром симметрии функции, она смещена по 
обеим координатам. Поэтому для малых интервалов определим значение одной промежуточной точки, для больших 
— двух точек. Результаты расчетов сведены в табл. 9. 
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Таблица 9 
Расчетные значения А„. в промежуточных точках 
У, км/ч Гц Кое У, км/ч Гц Кос 

70 8 0,1297 80 500 0,3188 
70 31,5 0,1597 80 2000 0,2586 
70 125 0,1968 80 8000 0,2098 
70 500 0,2429 87 500 0,2913 
75 8 0,1496 87 2000 0,2363 
75 31,5 0,1842 87 8000 0,1917 
75 125 0,227 93 500 0,2711 
75 500 0,2801 93 2000 0,2199 
80 8 0,171 93 8000 0,1784 
80 31,5 0,2105 100 500 0,2507 
80 125 0,2594 100 2000 0,2034 
80 500 0,3201 100 8000 0,165 
































На рис. 2 представлена графическая реализация полученных зависимостей. 
к 0,171 
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Рис. 2. Влияние скорости движения и частоты на величину коэффициента ослабления вибрации 


Выводы. Использование полученных зависимостей коэффициентов ослабления вибрации позволяет упро- 
стить расчет и повысить его точность октавных уровней вибрации и звукового давления структурной составляющей 
шума на рабочих местах локомотивных бригад при различных скоростях движения. Во-вторых, повышается его точ- 
НОСТЬ. 
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